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We consider strongly-correlated systems described by the multi-orbital Hubbard model in the 
atomic limit and obtain exact expressions for the chemical potential and thermopower. We show that 
these expressions reduce to the Heikes formula in the appropriate limits {ksT ^ U) and {ksT <^ U) 
and obtain the full temperature dependence in between these regimes. We also investigate the effect 
of a magnetic field introduced through a Zeeman term and observe that the thermopower of the 
multi-orbital Hubbard model displays spikes as a function of magnetic field at certain special values 
of the field. This effect might be observable in experiments for materials with a large magnetic 
coupling. 
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I. INTRODUCTION 

The thermopower of strongly correlated materials has 
received a lot of attention following the observation 
of a very large thermopower at room temperature in 
NaCo204 These and other similar materials can be 
described by single or multi-orbital Hubbard models with 
strong interactions. Several measurements of the ther- 
mopower of organic systems have been also carried out 
and the results were found to be describable by models of 
strongly interacting electrons on a lattice 0, ■ In these 
models, the hopping energies were assumed to be much 
smaller than the correlation energies and temperatures. 
This limit is the so-called atomic limit. The effects of 
strong correlations on the thermopower has been studied 
in the whole range described by the atomic limit earlier 
for the single band Hubbard model Q . The thermopower 
of the multi-orbital Hubbard model has also been studied 
earlier, but only in the low and high temperature limits of 
the atomic limit 0,0- In these limits, the thermopower 
is essentially the entropy carried per unit charge and can 
be calculated from thermodynamics alone without using 
the heat and charge currents of the system. A calculation 
over the entire atomic limit for the multi-band Hubbard 
model as done in this paper requires explicit forms for the 
charge and heat currents which are derived later. Other 
approaches to the calculation of the thermopower in these 
systems involve techniques such as dynamical mean-field 
theory 0, |1] and exact diagonalization |^ . 

Orbital degeneracy has a significant effect on the ther- 
mopower of a system. The increased entropy coming 
from the larger phase space due to the orbitals serves 
to enhance the thermopower. Such an effect has indeed 
been observed in LaCrOa and LaMnOs where the degen- 
eracy of the 3d orbitals causes an enhancement of ther- 
mopower 0, ^ . Orbital degeneracy has also been con- 
sidered as a candidate to explain the high thermopower of 
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the cobaltates [l^- The presence of many orbitals how- 
ever brings up issues of Hund's rule coupling and energy 
level splitting which complicate the study of thermoelec- 
tric transport. These effects are very difficult to take 
into account analytically if one is interested in calculat- 
ing the temperature dependence of the thermopower and 
have only been studied numerically 'l2]. The Hubbard 
model without Hund's rule coupling docs however lend 
itself to analytic treatment in the atomic limit as will be 
demonstrated. 

The effects of interaction on the thermopower beyond 
the atomic limit have been considered analytically for a 
single band Hubbard model in certain cases. Oguri and 
Maekawa studied the U ^ oo model near half fill- 
ing using a retraceable path approximation and Stafford 
[iJI used the Bethe ansatz to study the thermopower at 
low temperatures near half filling for the Hubbard chain. 
More recently, Paul and Kotliar Jl^ have calculated the 
thermopower in the Hubbard model in the infinite dimen- 
sional limit. The absence of an efficient large U perturba- 
tion theory for the Hubbard model makes the calculation 
of the thermopower and other transport coefficients dif- 
ficult in the presence of strong correlations. The atomic 
limit is thus most amenable to transport calculations in 
finite dimensions for a large range of temperatures and 
carrier densities. 

A magnetic field is expected to suppress the ther- 
mopower by lifting the degeneracy of levels and reducing 
the entropy. Such an effect has recently been observed 
experimentally in NaCo204 where the magnetic field 
lifts the spin degeneracy of the carriers drastically reduc- 
ing the thermopower. The effect is expected to be more 
pronounced in the presence of orbital degeneracy. The 
magnetic field affects not just the entropy but also that 
part of the thermopower which arises from a considera- 
tion of transport in the system. It will be shown that one 
sees spikes in the thermopower at certain special values 
of the magnetic field when that effect is considered. 

In section II, we treat the single band Hubbard model 
in a magnetic field and get exact expressions for the 
chemical potential and thermopower. We show that the 
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expression for the thermopower has particle-hole sym- 
metry and reduces to that obtained by a consideration of 
free spins in the U/UbT ^ 1 limit. Section III deals with 
the multiple orbital Hubbard model. Here we consider a 
system with degenerate orbitals and no Hund's rule cou- 
pling and derive expressions for the thermopower and 
fugacity as a function of temperature. We obtain the ap- 
propriate particle-hole symmetric Heikes expressions by 
considering the two limits U/ksT :s> 1 and U/ksT <^ 1. 
Finally, in section IV, we consider the multi-orbital Hub- 
bard model in the presence of a magnetic field. We notice 
the appearance of spikes for certain special values of the 
magnetic field and demonstrate how this follows as a con- 
sequence of energy conservation. 



Here g| = g and gi = —g. The thermopower is given by 

L2/Li+n/q 



where T is the temperature and fi the chemical potential. 
The coefficients L2 and Li are proportional to the Peltier 
coefficient and conductivity respectively and are given by 
the Kubo formulae. 

, _ /r drTr[e-/^(^-/-^)(Je Je(r) + J, J,(r))] 



II. SINGLE ORBITAL HUBBARD MODEL IN A 
MAGNETIC FIELD 



The Hamiltonian we wish to consider is 
H = -tJ24<^Cj+ia+h.c.+Uj2^nn3l+9Bj2i^:>^-^n) 

JO" j j 

(1) 

where t is the hopping parameter, U the interaction, g the 
Zeeman coupling and B, the magnetic field. The index 
J runs over all the sites in the chain. The Hamiltonian 
as written corresponds to a single orbital Hubbard chain 
with a magnetic field that couples to the carriers through 
only a Zeeman term. We note that all calculations in 
the atomic limit are carried out to lowest order in the 
hopping t (order t^) and to this order, the results are 
the same for higher dimensional lattices as well. We also 
note that for higher dimensional lattices, there is a term 
in the Hamiltonian arising from the magnetic field that 
couples to the hopping parameter. This term can also be 
neglected to lowest order in perturbation theory in the 
atomic limit. The expressions for the charge and heat 
current operators Je and Jq are given by 



i 

where n, and hj are the local charge and energy densities. 



a, the lattice constant and g, the charge of a carrier, 
the Hamiltonian given by Eqn. the currents are 
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For 
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L2 = 



Xr[e-/5(«-A'^)] 



(7) 



N is the number of carriers and (3 — l/fc^T. The time 
shifted current operators are given by 



(8) 



It should be noted that both Li and L2 are infinite in 
the atomic limit but their ratio is not. This is because 
the decays of the current-current correlation functions 
are assumed to be the same in both the heat and charge 
channels. This assumption presumably breaks down in 
the limit where real dissipation is introduced at higher 
orders of t/U. The procedure outlined from here on is 
similar to the calculation by Beni j3| in the absence of 
a magnetic field. We calculate the coefficients Li and 
L2 to lowest order in t (order t'^) and keep the hopping 
parameter in the calculation only in the current and not 
the exponential factors. We can also neglect the second 
term in the expression for the heat current Eqn. ^ The 
nearest neighbor hopping allows us to take the trace over 
only a pair of nearest neighbor sites (1 and 2) with the 
constraint of energy conservation n2-a = "i-o-- The 
resultant expression for the thermopower is 
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e'^^ is the fugacity and can be obtained from the 



equation 



dZ__ Z 
dx ^ X 



(10) 
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where Z = Tre^'^'^^''^^ is the partition function and p 
is the carrier density, with half-filling corresponding to 
p=l. In the present case, x is the root of a quadratic 
equation and is given by 
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{p - 1) cosh{f3gB) + V (1 - p)2 cosh^PgB) + p{2 - p)e-f^^ 
{2^p)e-PV 

I 



(11) 



The first two terms of the expression for the thermopower 
Eqnl^correspond to the contribution form transport and 
the third term corresponds to the entropy. The temper- 



ature dependence is given by the first two terms which 
disappear in the Heikes limit. In the limit of large ^r^, 
the thermopower reduces to 
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ha 

5 = — [Inp - ln(l - p) -ln{2cosh(/35B)} + l3gBt&n\i{(3gB)\ 
1 v ' v ' 

configurational spin 



for p < 1 and 



(12) 



S ^ — [ln(p - 1) - ln(2 - p)+ ln{2 cosh(/3.gB)} - PgB tanh(/3gB)] 

q " . ' ^ : — '■ ' 

configurational spin 



for p> I. 



(13) 



where the first term corresponds to the configurational 
entropy and the second term gives the spin entropy corre- 
sponding to free spins. The full temperature dependence 
is shown in Fig. 1. We can see that the thermopower 
is particle-hole symmetric and disappears at half filling 
{p = 1). For low and high values of the carrier density, 
the thermopower goes down in magnitude with increasing 
magnetic field (or decreasing temperature) corresponding 
to the reduction of entropy. For intermediate values, the 
thermopower changes sign going from particle like to hole 
like or vice versa because of the effect of the transport 
term. Indeed, the transport term can cause a change in 
the sign of the thermopower as a function of temperature 
even in the absence of a magnetic field . This effect is 
suppressed as one approaches the Heikes limit. 




III. MULTIPLE ORBITAL HUBBARD MODEL 

Let us now look at the multiple orbital Hubbard model. 
We assume that every site has a set of degenerate levels 
for the carriers and there is a probability for a carrier to 
hop from a level on one site to any level on an adjacent 
site. In addition, the carriers interact with all other car- 
riers on the same site with an interaction strength U . A 
more general model would include different interaction 
strengths for each set of orbitals and exchange energies 
(Hund's rule couplings). However such a model is not 
easy to treat analytically even at the atomic level and we 
consider this simplified model which still captures a lot 



FIG. 1: Thermopower as function of magnetic field for the 
single band Hubbard model. The curves are symmetric about 
half filling p = 1 and the thermopower goes down with applied 
field for low and high carrier concentrations. For intermediate 
values, the thermopower changes sign with field. 



of the essential physics. Our Hamiltonian is thus 



ja(3a ja 
+ ^ (14) 
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The currents are once again given by Eqn. 13 and in this 
case are 

"^e = 1^ ^ taf3Cj+laa'^jl3cr - C/3Cjq(jCj + 1/3<t (15) 



Ua 

Jq = -T 

in 



j a 0~f C7(T 



(16) 



r 



Here a, (3 and 7 are orbital indices while a and a' are where 
spin indices. The Kubo formulae Eqns. El and give the 
following upon reduction of the problem to a trace over 
two sites as before. 



L2 _ I2 



h = Tr[S({r^la<Tl,f^2/3<T2})2^ l^a/3| {n2aa -nif3a) 

aj3(T 



(17) and 



(18) 



I2 = Tr[B{{niaai,n2f3a2})U{^ \tal3\'^in2aa " nip„)'^n2a-a + ^ 1^7/3 P (^270-' " "■1/3<t' )^"2a<T}] 

and the traces have to be performed subject to the constraint that 



7//3rT' 



(19) 



(20) 



which is just the condition of energy conservation that 
says that the number of carriers on any two sites par- 
ticipating in a hop should be the same. We have used 
the fact that tap = t*^^ in the above expressions. In 
order to perform the traces, we make a further simplify- 
ing assumption that t^p = t for all a and f3. With this 
simplification, we are able to obtain the exact expression 



(21) 

where N is the number of orbitals per site and the fugac- 
ity X is given from Eqn. 1101 bv 



2N 



^ (T ) ^" " P)e-^^"'"-')/'a;" = (22) 

Tl=0 ^ ^ 

It is straightforward to show that the above equation has 
exactly one non-negative root for all values of p. For 



U ^ ksT, we obtain 



2iV(l-p) 



(23) 



which is positive only for p < 1 thereby preventing double 
occupancy. This gives us the multi-orbital Heikes formula 



S = — [ In p - ln(l - p) - hij^ - InjV] 

configurational ^pin orbital 



(24) 



It should be noted that Eqn. [21 is particle- hole symmet- 
ric about half filling p — N. For p < N, the transport 
part {L2/L1) goes to zero in the U/ksT 1 limit. For 
p > N, however, the contribution of the transport part 
has to be considered in addition to the entropic part and 
only a proper addition gives the right Heikes limit. For 
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FIG. 2: (Left) Thermopower plotted as a function of ksT /U 
for different values of p less than half filling for TV = 3. The 
thermopower goes to zero in the limit of small fcsT/f/ for 
p = 6/7 and is of opposite sign for values of p above and be- 
low this value. The curve for p — 1.2 diverges in this limit. 
(Right) Thermopower plotted as a function of p/2N for dif- 
ferent values of A'^. The thermopower is non-monotonic and 
oscillates with the number of maxima on either side of half- 
filling {p/2N = 1) equal to the number of orbitals. These 
oscillations go away in the limit of large ksT/U when all the 
curves collapse onto a single curve given by Eqn. 



p > N, we obtain 



S 



-^[ ln(2jV- 



p) -ln(p-2iV+l)- ln2 - IniV] 



configurational 



spin orbital 

(25) 

S is thus finite only for p < 2N — l\ip>N. For values 
of p between 1 and 2N — 1, the thermopower is divergent 
in this hmit except at half filling [p = N) where it is 
identically zero. In the opposite limit U <C fc^T, one 
obtains 



(26) 



2N - p 



The thermopower in this limit is 

5 = -^ln^ 

a p 



(27) 



where all values of p between and 2A'' give a finite ther- 
mopower. Eqn. is the so-called generalized Heikes 
formula previously obtained through direct combinatrics 



1^ . The expression for the thermopower in this limit 
is explicitly particle-hole symmetric since the contribu- 
tion to the thermopower for all values of p comes only 
from the entropic term. The temperature dependence 
between the two limits considered above is given by Eqn. 
I?T1 The thermopower changes sign as a function of fill- 
ing both below and above half-filling. Fig. 2 shows the 
dependence of thermopower on temperature and carrier 
concentration for = 3. Below half filling, the ther- 
mopower changes sign ai p = 2N/{2N + 1) which for 
A'^ = 3isp = 6/7, for small ksT/U and always remains 
negative for small ksT/U . Thus the curves with p > 6/7 
change sign as a function of the parameter kBT/U. It 
can also be seen that the curves are not monotonic at 
a fixed value of fc^T/f/ with p/2N and display oscilla- 
tions. The number of maxima on either side of half- filling 
{p/2N = 0.5) is equal to the number of orbitals. The os- 
cillations disappear in the limit of large kBT/U and all 
the curves collapse onto a smooth curve as given by Eqn. 



J 



IV. MULTIPLE ORBITAL HUBBARD MODEL is introduced through a Zeeman term, which couples dif- 
IN A MAGNETIC FIELD 

We now consider a multiple orbital Hubbard model in 
the presence of a magnetic field. Once again, the field 
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ferently to different orbitals. The Hamiltonian for such The charge and heat current operators are given by 
a system is 



jafla ja 
+ ^ '^jaa'rijpa' + YgaaC^aa'^]aa (28) 

^ I 



ih 



Jq 



Ua 



7^Q 



(30) 



r 



where the heat current now picks up an additional con- and 
tribution from the transport of magnetic energy in the 
same way as Eqn. 0] Once again, the thermopower is 
given by Eqns. ISl 1^ 171 and 1171 where 



h = Tr[B{{niaai,n2pa2}) Y |^a/3p(?^2a(T - "1/3<t)^ 



(31) 



h = Tl-[B{{niaai,n2f3a2}){Y \tap\^{n2aa ~ nifja)^ {Un2a-a- + QaaB) + ^ 1^7/31^(^270-' ~ nipa'Y n2aa}] 

after reducing the problem to a trace over two sites. The constraint on the trace is 

U{nip-a + ^ ni^„> - n2a-a - ^ n27<T') + {g[ja - gaa)B = (32) 

I 



This constraint differs from the one in Eqn. [5D| by 
the last term involving the magnetic field. This has in- 
teresting consequences as it produces spikes in the ther- 
mopower as shown in Fig. 4 at certain integer values of 
the magnetic field. This can be understood as follows: 
The constraint on the trace is just a statement of the 
conservation of energy. Without a magnetic field, all the 
orbitals are degenerate and the conservation of energy 
condition implies that the number of carriers on the two 
sites participating in a hop should be the same because 
the only energy involved is the Hubbard energy. In the 
presence of a magnetic field, the degeneracy of the lev- 
els is lifted and now the energy that has to be conserved 
is a combination of the Hubbard energy and the mag- 
netic energy. The constraint condition Eqn. 1321 reduces 
to Eqn. |2niexcept for special values of the magnetic field 
when {gprr — gacr)B /U is an integer. This involves two 



different kinds of orbitals {a and /3) where the Zeeman 
energy transferred is equal to the negative of the trans- 
ferred Hubbard energy. This transfer can take place only 
through two different kinds of orbitals because the hops 
are spin conserving and hence such spikes will not be ob- 
served in a system with just one kind of orbital. This is 
illustrated in Fig. 3. 

The spikes shown in Fig. 4 could be observed in exper- 
iments provided the magnitude of the Zeeman coupling is 
large. This is likely to happen in compounds with heavy 
transition and actinide metal ions such as the ruthenates 
and uranium based compounds where the strongly corre- 
lated electrons involved in transport are in the / orbital. 
The spikes in an actual experiment will be broadened due 
to thermal effects such as coupling with phonons or some 
other excitations which will relax the energy conserva- 
tion constraint. The effect of phonons on thermopower 
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in strongly correlated systems has been investigated ear- 
lier in a different context [T7|. 



4- 



> process that transfers Hubbard energy 



^ 

Energy transfer with multiple Ofbifals in a magnetic field 

FIG. 3: (Left) Energy conservation in the absence of a mag- 
netic field. The Hubbard energy has to be conserved and 
hence the total number of carriers on the sites participating 
in the hopping process has to be the same. (Right) Energy 
conservation with a magnetic field. The Hubbard energy and 
magnetic energy have to be conserved together and thus the 
number of carriers need not be the same on the two sites 
anymore. The difference in Hubbard energy can be compen- 
sated by the gain in magnetic energy because of the different 
coupling strengths of the orbitals involved in the hop. 



tained exact expressions for the thermopower and chem- 
ical potential and obtained the full temperature depen- 
dence of these quantities between the strong correlation 
U/ksT ^ 1 and high temperature U/kBT <^ 1 limits 
where they reduce to the corresponding Heikes formulae 
previously obtained using combinatric arguments. The 
thermopower of the multiple orbital system shows spikes 
as a function of magnetic field for certain special values 
of the magnetic field due to a transfer of commensurate 
Zeeman energy along with the Coulomb energy during 
hopping processes. We conjecture that it might be pos- 




FIG. 4: Thermopower as a function of magnetic field for the 
multiple band Hubbard model with A'^ = 3. The values of 
the coupling constant are gij = — 3/2,1721 = "1/2, 1731 = 
1/2, gi-^ = —1/2,(721 = 1/2,331 = 3/2. The thermopower 
consequently has spikes as a function of the magnetic field 
at special values of the field as explained in the text. The 
seemingly finite width of the spikes in the figure is due to the 
discrete nature of the points in the plot. 



sible to observe these spikes in some form in complexes 
of higher transition metal ions. 



V. CONCLUSIONS AND COMMENTS 
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